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Abstract. We present results for the equation of state of a graphene-like model in an effort to 
understand the properties of its quantum phase transition. The Nf fermion species interact 
through a three dimensional instantaneous Coulomb potential. Since there are no reliable 
analytical tools that work for all values of Nf and the coupling constant g , we rely on Monte 
Carlo simulations to calculate the critical properties of the model near the phase transition. 
We consider the four-component formulation for the fermion elds, which arises naturally as 
the continuum limit of the staggered fermion construction in (2+1) dimensions. In the limit 
of infinitely strong Coulomb interaction, the system undergoes a quantum phase transition at 
a critical number of fermion species Nfc ~ 4.7. We also calculate the values of the critical 
exponents at the quantum phase transition. 



1. Introduction 

There has been considerable recent interest in graphene sparked by its discovery and subsequent 
experimental study. In brief, for a carbon monolayer having one mobile electron per atom, a 
simple tight-binding model shows that the spectrum of low-energy excitations exhibits a linear 
dispersion relation centred on zeroes located at the six corners of the first Brillouin zone (see 
e.g. [T]). Using a linear transformation among the fields at two independent zeroes it is possible 
to recast the Hamiltonian in Dirac form with Nf = 2 flavors of four-component spinor tp, the 
counting of the massless degrees of freedom coming from 2 carbon atoms per unit cell x 2 zeroes 
per zone x 2 physical spin components per electron. Electron propagation in the graphene layer 
is thus relativistic, albeit at a speed vp ~ c/100. Although this is a factor of 100 slower than 
the speed of light in vacuum, it is still much faster than the speed of electrons in an ordinary 
conductor. The implications for the high mobility of the resulting charge carriers (which may be 
negatively-charged "particles" or positively-charged "holes" depending on doping) is the source 
of the current excitement. The stability of the zero-energy points is topological in origin, as 
emphasised by Creutz [2]. 

While the above considerations apply quite generally, a realistic model of graphene must 
incorporate interactions between charge carriers. One such model, due to Son [3j, has Nj 
massless fermion flavors propagating in the plane, but interacting via an instantaneous 3d 
Coulomb interaction. In Euclidean metric and static gauge d^Ao = the action reads 





where e is the electron charge, V = Aq, and the 4x4 Dirac matrices satisfy {7/^,7^} = 25^^, 
fj, = 0,1, 2, 3. In our notation x is a vector in the 2d plane while the index i runs over all three 
spatial directions. In the l&ige-Nf limit the dominant quantum correction Il{p) comes from a 
vacuum polarisation fermion- antifermion loop and the resummed V propagator becomes 



A(,) = (o„-'(p)-n(p))-= f|! + iiiJ^ , (2) 




where = {po,p)'^ = Po + Vp\p\'^ , and Dq{p) corresponds to the classical propagation of V. 
In either the strong coupling or large- A'^y limits Di{p) is thus dominated by n(p), the relative 
importance of the original Coulomb interaction being governed by a parameter A = |n/Z?o|po=o- 

The chiral symmetry breaking, due to the spontaneous condensation of particle - hole pairs, 
is signalled by an order parameter (ipip) ^ 0. Physically the most important outcome is the 
generation of a gap in the fermion spectrum, implying the model is an insulator. Son [3] 
postulates that this insulating phase exists in the corner of the phase diagram corresponding to 
large and small A'^^, and in particular that the insulator-conductor phase transition taking 
place at Nf = Nfc in the strong-coupling limit — > 00 is a novel quantum critical point. The 
value of Nfc, and the issue of whether it is greater or less than the physical value Nf = 2, must 
be settled by a non-perturbative calculation. A recent estimate obtained by a renormalisation 
group treatment of radiatively-induced four- fermion contact interactions, is Nfc = 2.03 [4J. 

The proposed physics is very reminiscent of the three dimensional Thirring model, which 
is analytically tractable at large Nf, but exhibits spontaneous chiral symmetry breaking at 
small Nf and strong coupling [Sj [6]. Arguably the Thirring model is the simplest field 
theory of fermions requiring a computational solution: the location of the phase transition 
at Nf = Nfc in the strong coupling limit has recently been determined by lattice simulations 
to be Nfc = 6.6(1) [7]. The apparent similarity of the two systems has led us to propose a 
Thirring-like model pertinent to graphene, with Lagrangian 



S2 = ^ / dxo(fx 

a=l 



(3) 



As for ([T]) we assume a large- A^j limit to estimate the dominant vacuum polarisation correction; 
the resultant propagator for V is 




D2{P)= -, + ^-^ . (4) 



In the strong-coupling or large- A'^y limits, D2 coincides with Di, implying that the fermion 
interactions are equivalent. It is also the case that liuip^oo D2{p) = limx-^oc Di{p). This last 
limit is important because critical behaviour in the Thirring model is governed by a UV-stable 
fixed point of the renormalisation group We thus expect predictions made with the model 
([3]), and in particular critical behaviour such as the value of Nfc, to be generally valid for Son's 
model ([1]) in the limit of large A. 



2. Lattice Simulations 

The lattice model studied in this paper is related to the lattice Thirring model studied in 
[5], except that in this case we used the compact formulation for the auxiliary boson field. 
The simulations were performed with the standard hybrid molecular dynamics algorithm. The 
novelty of [7J is that it is the first study of the Thirring model by lattice simulations in the 
strong-coupling limit ^ 00. Since we aim to repeat the strategy here we should discuss how 



0.2 



0.1 



j-i D Q O D 



16^m=0.01 
° 24\j=0.01 
° 24^m=0.02 

32^ m=0.01 



Figure 1. (V'V') vs. 1/^^ near ^peak for Nf = 2. 



this was done. First, note that the vacuum polarisation calculation leading to the results 1^ 
does not go through in quite the same way for the lattice regularised model; rather, there is an 
additive correction which is momentum independent and UV-divergent: 



n'-«(p)=n-'^^(p)+gV(m), (5) 

where J(m) comes from incomplete cancellation of a lattice tadpole diagram [5j. This extra 
divergence not present in the continuum treatments can be absorbed by a wavefunction 
renormalisation of V and a coupling constant renormalisation 

In the large-A^/ limit we thus expect to find the strong coupling limit of the lattice model at 
g% ^ oo implying g^^. For g^ > g^^ £'iatt(p) becomes negative, and 5iatt no longer 

describes a unitary theory. 

Away from the large- A^/ limit, where chiral symmetry may be spontaneously broken, there 
is no analytical criterion for identifying gfi^j^; however in this case a numerical calculation of 
{^l'■^p) shows a clear peak at g^ = 5peak' whose location is approximately independent of both 
volume and m, indicating an origin at the UV scale [7j. Fig. [T] exemplifies this behaviour in 
the lattice model with Nj = 2 on system volumes L^: for L > 24 we identify l/^p^ak — 0-3- 
Since for orthodox chiral symmetry breaking the magnitude of the condensate is expected to 
increase monotonically with the coupling strength, we interpret the peak as the point where 
unitarity violation sets in, ie. gf^^ ^ Speak- simulations performed at g^ = gf^^ to 

explore the strong coupling limit, and find clear evidence for a chiral symmetry restoring phase 
transition at a well-defined iV/c- As shown in Fig. [2] l/f^p^^j^j^ decreases as Nf increases from 2 to 
4.75. Near Nf ~ 4.75 the curve reaches a minimum, implying a significant change in the strong 
coupling behavior of the model. The rapid change in the value of the condensate at Nf ~ 4.9 
implies a signifacant change in the strong coupling behavior of the model. This can be a phase 
transition that separates a chirally broken from a chirally symmetric phase. Since our model has 
anisotropic interactions in the spatial and temporal directions the correlation lengths near the 
transition diverge with different exponents, i^s in the spatial directions, and vt in the temporal 
direction. In this case a modified hyperscaling relation holds: 



ut + {d- = 7 + 2/3. 



(7) 
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Figure 2. l/gl^^^ vs. Nf. 



Figure 3. Equation of State fits to data 
generated witli different fermion masses on 
different lattice sizes. 



We carefully monitored the finite size effects in the spatial directions by comparing the values 
of the condensate at m = 0.01, fl'peak(-^/) extracted from simulations on 16^ X 48 and 24^ x 48 
lattices. Our results show that finite Lg effects are negligible for Lg > 16. We fitted the data 
generated on several lattices with Lt = 16, Lg = 48,64 and m = 0.01, ...,0.04 to the following 
RG-inspired equation of state 

m = A[{Nf - Nf,) + CL~^']{i>i,Y + B{i^^p)', (8) 

where p = 6 — (i^^ ■ The results show clearly that the model undergoes a second order quantum 
phase transition and the values of the most significant parameters are: Nf, = 4.74, 6 = 3.55, 
p = 822, and ut = 2.56. The data and the fitted curves are shown in Fig. [3l The value of 
Nfc = 4.74 is close to the value of Nf ^ 4.9, where a rapid change in the l/Speak occurs. 

3. Summary 

We presented results from a non-perturbative Monte Carlo study of a Thirring-like model 
pertinent to graphene. Our analysis of its equation of state at strong coupling reveals that 
the model undergoes a second order quantum phase transition at Nfc ^ 4.74, which is greater 
than the physical Nf = 2 value. This result implies that freely suspended graphene may be an 
insulator for couplings g larger than a critical coupling g,- 
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